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The statistical mechanics of a system of non-relativistic 
charged particles in a constant magnetic field is discussed. 
The spatial dimension D is arbitrary with D > 3 assumed. 
Calculations are presented from first principles using the ef- 
fective action method. For D > 5 the system has a phase 
transition with a Bose condensate. We show how the effec- 
tive action method method deals with in a very natural way 
with the condensate and study it's role in the magnetization 
of the gas. For large values of the magnetic field we show 
how the magnetized gas in D spatial dimensions behaves like 
the free Bose gas in {D — 2) spatial dimensions. Even though 
for D = 3 the magnetized gas does not have a phase tran- 
sition for any non-zero value of the magnetic field, we show 
how the specific heat starts to resemble the result for the free 
gas as the magnetic field is reduced. A number of analytical 
approximations for the magnetization and specific heat are 
given and compared with numerical results. In this way we 
are able to study in precise detail how the B ^ limit of the 
magnetized gas is achieved. 



I. INTRODUCTION 

Bose-Einstein condensation (BEG) has been an 
integral part of physics for many years. From a theo- 
retical standpoint this phenomenon has been intensively 
studied as a possible means of explaining the superfluid 
transition in He* j|] , and also as an experimentally real- 
isable possibility in dilute atomic gases. Recently BEG 
has been achieved in atomic systems of dilute alkali metal 
atoms and for this reason a great deal of interest 

is now being focussed on these systems. 

Another system worthy of study is the condensation 
(or lack thereof) of charged bosons in a homogeneous 
magnetic field. Initial work on the nonrelativistic case 
1^ was substantially improved on by Schafroth and 
was then generalized to dimensions other than three by 
May The primary motivation for these efforts 

was to formulate a plausible theory of superconductivity. 
Although this attempt failed as far as standard low Tc 
superconductors are concerned, interest in the subject 
has revived as explanations of high Tc materials are now 
required [ p^ . 

More recently, several detailed studies have been made 
of the ideal system of nonrelativistic charged bosons, 
both numerically [|l3| , and from a theoretical point of 
view in both two [HTand three lEsNlql dimensions. Also 



a detailed study of the boson gas in three dimensions with 
trapping harmonic potentials and the presence of a mag- 
netic field has recently been made using the path-integral 
formalism jl^ . These authors find comparable behavior 
to that found using Mellin transforms, and a comparison 
between the two methods is very instructive. However 
due to constraints of space, the contrast between the two 
methods will not be discussed here in depth. 

The behavior of the relativistic case has also been dis- 
cussed In addition more realistic systems 
which may occur in high Tc superconductors |23|^25|| have 
also been examined. 

The outline of our paper is as follows. In section || we 
give a brief description of the free Bose gas in D spatial 
dimensions. We will later compare the analogous results 
for the magnetized gas to these free-field ones. Section 
III presents the effective action method and and applies 
it to the magnetized gas. We concentrate initially on the 
specific heat and show how the presence of the magnetic 
field alters the behavior from that found for the free gas. 
We study numerically what happens for large and small 
magnetic fields. Results are also obtained for the mag- 
netization. In Section IV we describe how it is possible 



to obtain analytical results for the critical temperature 
(when -D > 5), the magnetization, and the specific heat 
when the magnetic field is weak. We are able to analyt- 
ically confirm the numerical results concerning the limit 
_B — !■ 0. For the D — 3 gas we give approximations valid 
at the critical temperature for the free Bose gas. (Pre- 
vious approximations were only valid for temperatures 
larger than this.) A number of appendices contain de- 
tails of the expansions used for the analytical results. 



II. BOSE-EINSTEIN CONDENSATION WITH NO 
EXTERNAL FIELD 

In this section we wish to review very briefly some of 
the basic properties of the free Bose gas. The spatial 
dimension D will be arbitrary. For D = 3 the analysis 
is standard textbook material p6|-p9[|. The absence of 
Bose-Einstein condensation when 13 = 2 is also widely 
known |3C|,[l0[|. Spatial dimensions D > 3 have also been 
studied [^,^. 

We will consider an ideal gas of N spinless bosons con- 
fined in a large box of volume V in D spatial dimensions. 
The infinite volume limit will be taken with Af = ^ fixed 
as is conventionally done. The energy levels of the system 
are 
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(2.1) 



before the infinite volume limit is taken if we impose 
periodic boundary conditio ns o n the field. We set h — 1 
throughout the paper. In (2.1), rii = 0, ±1,±2,... and 
Li is the length of the box in the ith direction. In section 
[II we will show how the effective action formalism can 
be used to study this problem (see references p6| , p3| , p4[ 
for reviews). For now we stick with the conventional 
thermodynamic expressions. 
The internal energy is given by 
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(2.2) 



where fji is the chemical potential and (3 = in units 
with the Boltzmann constant fc^ = 1. The grand canon- 
ical ensemble is used here. 

By taki ng t he infinite volume limit, we can replace the 
sums in (2.2) over rii with integrals. After expanding 
the denominator of ( ^.2[ ) in a geometric series, and using 
(2.1) for the energy levels, the integrals over rii may be 
performed with the result 
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We have defined the polylogarithm function Li^ [z] by 
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We note the property (assuming p > 1) 
Lip [1] = Ch(p), 



(2.4) 



(2.5) 



where Ci, (p) denotes the Riemann ^-function. 

The essential feature of Bose-Einstein condensation as 
a phase transition is exhibited in the behavior of the spe- 
cific heat. The specific heat at constant volume Cy is 
defined by 
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The particle number N, given by 
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(2.6) 



(2.7) 



is held fixed when computing the derivative in (2.6). The 
chemical potential is not fixed. From (2.2), (|2.6D we find 
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Here ' denotes the derivative with respect to (3 holding 
V and N fix ed. To calculate (/3/i)' we differentiate both 
sides of (2.7) to give 
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(2.9) 



Substitution of (2.9) into ( |2.S| ) and use of ( p.7| ) shows 
that 
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If we confine ourselves to D > 3, as we do for the rest 
of the paper, then a critical temperature Tq can exist at 
which the chemical potential /i = 0. From ( ^.7|) we find 
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For T < Tq the chemi cal po tential remains fixed at /i = 0. 
It is easy to see that ( ^.1C| ) only holds for T > Tq. When 



T < Tq, the specific heat may be evaluated from (2.8) by 



setting the term with (Pfi)' to zero and also fi to zero. 
This results in 
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for T < Tq. 

The two expressions (2.1C),(2.12) may be used to com- 
pute the specific heat for all temperatures, and also to 
study the behavior at T = T q. It is easy to see from the 
polylogarithm function (2^) that if p < 1, Li p \z] — > oo 
as z ^ 1. This means that the second term of ( 2.10| ) van- 
ishes as /i ^ for D = 3, 4, but is finite for D > 5. Thus 
the specific heat is continuous at T = Tq for D = 3, 4, but 
discontinuous at T = Tq for D > 5. The discontinuity 
for Z? > 5 is easily computed in terms of Riemann l^- 
functions. The continuity of Cy for D = 3 is well known 
iel-Egl. The behavior for D > 3 can be found in il|,0. 



III. BOSE-EINSTEIN CONDENSATION IN A 
CONSTANT EXTERNAL MAGNETIC FIELD 

When a magnetic field is applied to a gas of charged 
bosons in three spatial dimensions the energy spectrum 
(in the infinite volume limit) contains a discrete harmonic 
oscillator-like part as well as a continuous part. The dis- 
crete part is just the Landau level quantization For 
D > 3 there may be a number of discrete components 
because the magnetic field is not described by a vector, 
but by an antisymmetric tensor with more than one in- 
dependent component ]3^ ] . For simplicity we will restrict 
our attention to the case of only a single nonzero compo- 
nent in the present paper. We wish to provide a treat- 
ment similar to that for the free Bose gas when a nonzero 
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magnetic field is present. In particular we will study the 
specific heat and see how the presence of a magnetic field 
alters the behavior from that found for the free Bose gas 
in section ||. Also we will compute the magnetization 
and study the Meissner-Ochsenfeld effect in detail. The 
formalism used is the effective action method as reviewed 
in references [^^. This formalism allows the nonzero 
condensate (if there is one) to be treated in a very natural 
manner. 



A. Thermodynamic potential and phase transitions: 
General formalism 

The thermodynamic potential is usually defined by 
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However, in the effective action method there is another 
term present if there is a nonzero condensate described 
by a background field ^ . This is 
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Here = V^* — ieA^ is the usual gauge-covariant 
derivative. The complete thermodynamic potential is 



T=iO- 



(3.3) 



(Actually, if we are interested in the dynamics of the 
magnetic field there will be an additional term involving 
a Maxwell action. We will consider this in section [HD 



below). Given the thermodynamic potential, all quanti- 
ties of interest can be calculated. 

The presence of a condensate is signalled by a nonzero 
value for ^ . This is associated with symmetry breaking 
as discussed in the relativistic case • In our case ^ 

must satisfy 
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We can solve this by expanding \E'(a;) in terms of the 
stationary state solutions to the Schrddinger equation: 
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(3.5) 



If we write 



*(a;)===^C7„/„(a;) (3.6) 

n 

for some coefficients C„, and assume that the set of so- 



lutions fn{x) forms a complete set, then (3.4) results in 
= (S„ - etJL)Cn. (3.7) 
We will define a critical value of /i, say ^c, by 
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(3.8) 



where Eq is the lowest energy level. If /i < /^c, then the 
only solution to ( p/7| ) is for C„ = 0, which corresponds 
to 4" = 0. There is no condensate in this case associated 
with symmetry breaking and a phase transition. However 
if fj, can reach the value Hc defined in (3.8) for some 
temperature T^, then Cq in (3.7) is undetermined and 



we can have a nonzero condensate described by 
^{x) = Cofoix). 



(3.9) 



The temperature Tc at which fi ~ '^s called the critical 
temperature. 

For the case of the free gas considered in section |l| we 
have Eq — 0, so that /ic — 0. The critical temperature 
Tc is then the value of the temperature at which the 
chemical potential vanishes as stated earlier. If the spa- 
tial dimension 13 > 3 a critical temperature exists and 
signals a phase transition with a nonzero value for (5* 
is constant for the free Bose gas). Associated with this 
phase transition is a growth in the number of particles in 
the ground state. 

For some systems it is possible to have a sudden growth 
in the occupancy of the ground state without a phase 
transition. In this case fj, never reaches the critical value 
of fic but instead approaches it asymptotically. The 
speed at which the ground state particle number builds 
up depends on how fast fi approaches fic- Because fi 
never reaches we have ^E" = 0, and no symmetry break- 
ing. As we will see below, for a charged Bose gas in a 
constant magnetic field the spatial dimension D deter- 
mines whether or not a phase transition occurs. 

We can use our expression (3.3) for r2 to find the total 
charge Q, since Q = — with V, B,P, and ^ held fixed. 
It is convenient to write 



Q = Qo + Qi 



where 
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if we use ( |3.2D , ( |3.9| ), and 
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From (O) we find 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



If we can always solve Q — Qi for /i for all temperatures, 
then ^' = 0. There is no condensate, symmetry breaking, 
or phase transition in this case. If it is not possible to 
solve Q = Qi for /i, then we must have Qo and find 
a nonzero value for 
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B. Thermodynamic potential: Constant magnetic 
field 



The formalism outlined in section III A will now be 



applied to the £)-dimensional charged Bose gas in a con- 
stant one-component magnetic field. We will assume 
D > 3 here and pick the mag netic field in the z-direction. 
It is possible to solve (3.5) for the energy levels and 
corresponding eigenfonctions |3^]. We have (choosing 



where n = 0, 1, . . . labels the Landau level, and ki = 
0, ±1, ... if we impose periodic boundary conditions on a 
box as in section ^ We have defined 

Lu^—. (3.15) 



The energy levels (3.14) are degenerate with degeneracy 
eBLiL2 
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As in section ^ we will be interested in the large box 
limit with — > cxd. In this limit we can replace the 
sums over ki resulting when (3.14) is used in (3.1) with 
integrals. A change of variables gives 
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This may be evaluated by expanding the logarithm in 
its Taylor series and then performing the integral over k. 
We find 
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where the second line has followed by performing the sum 
over n. 

At this stage it is useful to define a dimensionless tem- 
perature. We define 



(3.19) 



X is seen to be the ratio between the energy gap between 
successive energy levels uj and the thermal energy k^T. 
The lowest energy level from (3.14) is i?o,o = 7- From 
( |3.8|) the critical value for fi is 



(3.20) 



We will define a dimensionless chemical potential e by 



e/i = [ i - e 



(3.21) 



A phase transition is characterized by a cr itical temper- 
ature Tc at which e = 0. The expression ( |3.18 ) may be 
written in terms of the dimensionless variables x and e: 
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It is convenient to introduce some compact notation for 
the class of sums we will encounter in order to simplify 
formulae. Let 
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(3.23) 



With this notation we may write ( p.22[ ) as 

^T^o = -ioVi—j ^,hD,0]. (3.24) 



Various thermodynamic quantities involve derivatives of 
the thermodynamic potential. We will initially consider 
the charge. 



From (3.12), using (3.24) we have 

d 



(3.25) 



Whether or not a phase transition occurs is determined 
by the convergence or divergence of Si [2 — D,0] as e ^ 
0. If this sum diverges as e — > then we will always 
be ab le to solve Q = Qi for /i. As discussed in section 
III A this means that there is no phase transition. From 
(|3.23 ) it is easy to see that [a, 0] diverges as £ ^ 
for a > —2. With a — 2 — D, this means that there is 
no phase transition for Z) = 3, 4 dimensions. For D > 5, 
Si [2 — D, 0] converges as e ^ and there is a phase 
transition. In this case {i.e. for D > 5) there is a nonzero 
condensate characterized b y ^ 7^ 0. 

For 13 > 5 we may use (3J) to find 4". If we choose 
the gauge 



Ai = -By, A2 = ... = Ao=0, 
for the vector potential, then 



(3.26) 



(3.27) 



is the eigenfunction corresponding to the lowest energy 
Eq — ^. a is a c ons tant ch osen to normalize /q. 1^'"' 
is found by using (3.9) with (3.27) in the general expres- 
sion ( |3.2| ). We will retu rn to this when we discuss the 
magnetization in section [II D| . 
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C. The specific lieat 



The internal energy is given by U 



and held fixed. Using (3.1 — 3.3) we find 
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{(30.) with 



The first term accounts for any nonzero condensate, and 
the second term is easily seen to be 
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which is t he u sual expression for the internal energy (see 
equa tion ( ^.2|) with a different definition for /x). With 
(3.5) and ( |3.9D we find 



U^Eo\Co\^ + -^mT^o) 



(3.30) 



If Eq = 0, then the contribution from fJ'^*'^ to the internal 
energy vanishes, and the i nternal energy is given by the 
standard expression (3.29). This is the situation for the 
free Bose gas discussed in section ||. For the constant 
magnetic field, Eq = ^, so that if 7^ we must in- 
clude the condensate contribution to obtain the correct 
expression f or the energy. 



If we use ( 3.11 ), then ( 3.30 ) may be written as 
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With (3.24) for flr^o, and differentiating with respect to 
(3 keeping lu, V fixed results in 
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Noting (3.25) allows us to write this as 
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Subst itution of ( 3.32 ) back into ( 3.31 ) and noting from 
(3.10) that Q = Qo + Qi where Q is the total charge, 
results in 
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The ex press ion for the internal energy we have just ob- 
tained ( |3.34 ) holds whether there is a nonzero condensate 
(^ 7^ 0) or not. In cases where a phase transition does 
occur, the role of the condensate is crucial for obtaining 
the correct expression for the internal energy. If w e had 
neglected the contribution coming from then (H^H) 
would have Qi in place of Q. This would then lead to 
an erroneous expression for the specific heat, since Q is 
fixed whereas Qi is not. 

The specific heat at constant volume was defined in 
( p!6| ). The quantities held fixed are Q, and w = 
when the differentiation is performed. This means that 
the first term in (3.34) makes no contribution to C^. Just 
like the free Bose gas discussed in section ||, we must 
distinguish between the expressions above and below the 
the critical temperature if there is a phase transition. For 
D > 5, the critical temperature Tc is defined by 
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Here Xc = Pc^ with Pc = T~^, and the sum E^ [2 — £), 0] 
is evaluated with e = and x — Xc- Unlike for the free 
Bose gas, it is not possible to solve for Tc analytically. 
We will _j^turn to an approxi mate evaluation of Tc in 
We have solved ( 3.35 ) numerically to find 



IV 



section 

Tc- The results for small values of uj are shown in figure 
|l|. As expected, when w ^ we have Tc —>Tq. 




FIG. 1. For D > 5 this shows the value of /3c compared 
to the free Bose gas critical value of Po. For D = 3, 4 be- 
cause there is no critical temperature when a magnetic field 
is present we have plotted the inverse temperature of the spe- 
cific heat maximum in units of Pq. 

For small values of cj, Pc is close to /3o (an approximate 
analytical expression will be given in section IV). As the 
strength of the magnetic field is increased the deviation 
of Tc from Tq becomes more pronounced, but lessens 
with increasing dimension. In all cases where a phase 
transition occurs, the critical temperature is lower than 
the result for the free Bose gas as observed in |]39| . 

Let denote the specific heat fo r T > Tc and be 
the specific heat for T < Tc- Using ( 3.34| ) we have 
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(3.36) 



in a constant magnetic field is continuous at the critical 
temperature for £' = 5,6 and discontinuous for D >7. 

In the cases Z? = 3,4 where there is no phase transi- 
tion, e never vanishes. The specific heat is given by ( 3.3j ) 
in these two cases for all temperatures. When £' = 3,4 
the specific heat is a perfectly smooth function of tem- 
perature. 

Graphs showing the specific heat for D = 3 — 7 are 
shown in figures — pi 



The terms in (^^^§^^ arise from the fact that the charge 

Q rather than the chem ical potential /i is held fixed. By 
setting Q = Qi in ( 3.25| ) (valid since T > Tc), and differ- 
entiating both sides with respect to (3 with Q held fixed 
gives 



fix 



d{ex) 
dp 



Si [4 - D, 0] 
D-2 



xT., [4 - D, 1] 
J:,[2-D,0]). (3.37) 



Substitution of ( 3.37 ) into ( 3.36| ) leads to 



C> = xV i^—j i^2x^^, [4 -D,2]+ x^S, [4 - D, 1] 

+ {D - 2)xY.^ [2 -£>,!] + ^^^^ ^'' Si [-£1, 0] 
x{D - 2)S2 [4 - £>, [2 - £), 0] 
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xH^, [4-i?,l])2 
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(3.38) 



For T < Tc we have e = fixed. (Equivalently, M =^ 
is fixed). The result of this is that is given by (3.3(;) 
but with the last two terms set to zero, and £ = in all 
other terms: 



+x^T,^ [4 - £), 1] + (£1 - 2)j;E2 [2 - £>, 1] 
£)(£)- 2) 



-xEi [-£',0] 



(3.39) 



£ = 



By comparing in ( ^.38 ) with C^f in ( 3.39 ) it can be 
seen that whether or not the specific heat is continuous 
at the critical temperature is determined by the behavior 
of El [4 — 13, 0] as e ^ 0. The two expressions will only 
agree if E, [4 ~ D, 0] diverges in this limit. From the 
definition ( 3.23 ) this only happens for £' = 5,6 (Recall 
that we are assuming D > 5 here so that Tc exists). We 
conclude that the specific heat for the charged Bose gas 
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FIG. 2. This shows the specific heat at constant volume 
in units of the total number of particles as a function of the 
inverse te mpe rature in units of /3o where /3o = Tg'^ with To 
given in ( 2.11 ) for D — 3. For comparison the free gas result 
is also shown {u — 0). 
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FIG. 3. The specific heat at constant volume in units of the 
total number of particles for the magnetized gas with £> = 4. 
The top curve (cj = 0) is the result for the free Bose gas. 
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FIG. 4. The specific heat at constant volume in units of the 
total number of particles for the magnetized gas with D — 5. 
The curves are given as a function of /? = in units of /3o 
rather than /3c- The discontinuous free gas result is labeled 




FIG. 5. The specific heat at constant volume in units of the 
total number of particles for the magnetized gas with D — 6. 
The discontinuous free gas result is shown as tj = 0. 




FIG. 6. The specific heat at constant volume in units of the 
total number of particles for the magnetized gas with D — 7. 
This is the first spatial dimension for which the specific heat 
has a discontinuity at the critical temperature. The free gas 
result [lj — 0) is also shown . 



For the = 3, 4 gases, the presence of the magnetic 
field is seen in figures to round olT the familiar sharp 
behavior at the transition temperature for the free gas. 
In both cases as the magnetic field is reduced the curves 
for the specific heat tend towards the free-field results. 
The maximum in the specific heat becomes sharper, and 
the temperature at which the maximum occurs tends to- 
wards the value Tq as the magnetic field is reduced. For 
B Q, although there is no phase transition character- 
ized by a critical temperature and a nonzero condensate, 
as B is reduced the specific heat starts to look more and 
more like the free gas result. We will look at this an- 
alytically in section IV. In both cases, as /3 ^ (or 
T oo) the specific heats for D — 3,4: approach the 
classical Maxwell-Boltzmann results of 1.5 (for D = 3) 
and 2 (for D = 4). 

The behavior of the specific heat as B is increased is 
also of interest. Increasing the value of B tends to reduce 
the specific heat maximum and to broaden the curves. In 
fact for large values of B if we examine figure it can 



be seen that the curves approach the value of^ as B 
is reduced before rising sharply to the classical result of 
|. This demonstrates that the specific heat for the gas 
with D = 3 resembles the specific heat for the free gas 
with Z? = 1 in a strong magnetic field. The value of i 
is the classical value for the one-dimensional gas. For 
the magnetized gas with D = 4, figure ^ shows that the 
specific heat approaches the value of 1 in a strong mag- 
netic field before the sharp rise to the classical value of 
2. Again this shows the reduction in the dimension by 



1 



^ . as /3 ^ for the free Bose gas in two 
spatial dimensions. This is totally consistent with the 
approach used in references |^,^ in which the leading 
behavior of thermodynamic quantities was studied in a 
general setting by using the lowest energy solutions. As 
B is increased, the gap between the ground state and the 
excited states becomes larger; thus the leading contribu- 
tion would be expected to come from the ground state. 

The results for the specific heat of the 5-dimensional 
gas are shown in figure ^. In this case the specific heat 
for the free Bose gas is discontinuous. Nevertheless as 
the magnetic field is reduced the specific heat curves 
approach the free gas result. The peaks of the specific 
heat start to become sharper and the slope of the curve 
steeper. The classical Maxwell-Boltzmann result of | is 
reached as /9 ^ 0. Just as for the cases Z? = 3, 4, the gas 
exhibits a reduction in the effective dimension for large 
values of the magnetic field. This time we would expect 
to find the specific heat curves looking more and more 
like the familiar form for the specific heat of the free gas 
in three spatial dimensions. For clarity we have given the 
large magnetic field results on a separate graph in figure 
1^ over a large temperature range. 

The reduction in dimension should be evident with the 
specific heat approaching the value |, the classical free 
gas result for = 3, before rising to the value | as 
f3^0. 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 

/S//3. 

FIG. 7. The specific heat is shown in units of for the 
magnetized gas with D = 5 for large values of the magnetic 
field. The curves resemble the free gas specific heat in three 
spatial dimensions over a wide temperature range (compare 
with tj = in figure |^ . 

The specific heats for = 6, 7 are shown in figures |^,^ 
respectively. For D = 6, we see the peaks sharpen as the 
free gas result is approached. Like the 5-dimensional gas 
the curves are always continuous, but the slope becomes 
greater as B is reduced. For strong magnetic fields the 
apparent reduction in dimension again takes place with 
the specific heat resembling the free Bose gas in four spa- 
tial dimensions. The results for D = 7 are shown in fig- 
ure ^, and are included because this is the first spatial 
dimension for which the specific heat in the presence of 
a magnetic field is discontinuous. For strong magnetic 
fields the specific heat resembles the free 5-dimensional 
gas, and for weak fields the free Bose gas limit is ap- 
proached. 



D. Magnetization 

Even though for D = 3 there is no phase transition 
which can be associated with Bose-Einstein condensa- 
tion, by studying the magnetization of the charged gas 
Schafroth showed that the gas exhibited the Meissner- 
Ochsenfeld effect. The generalization to other spatial di- 
mensions was performed later [ |10yil|jr^ - p8| , ^. W e will 
show how the formalism described in section III A can be 



used to obtain the magnetization. In particular, the role 
of the condensate 5* for > 5 will be examined carefully. 

The simplest way to see the effects of magnetization 
is by studying how the field equations for electromag- 
nctism arc affected. To do this we must include a term 



in the thermodynamic potential (3.3) for the electromag- 
netic field. We will use Heaviside-Lorentz rationalized 
units as usual in quantum field theory. (A discussion of 
the various units and how this alters the expression for 
the magnetization is given in | |4l[ |. Of course the physics 
of the situation should be independent of this arbitrary 
choice). We add on 



: p. . p«j _ 7' 4 



(3.40) 



to (pj 



where Fij — diAj 



is the field strength 



tensor describing the magnetic field, and Jl^^ is the ex- 
ternally applied current which is responsible for setting 
up the magnetic field. The complete thermodynamic po- 
tential is 



T#0 



(3.41) 



where 17^°) and r^T/o are given in (3.1) and (pT^). 

Variation of f2 with respect to the magnetic field 
results in 



P) . TJV — p 



where 



jji] ^ pi] 



SF; 



(3.42) 



(3.43) 



if*-' is the Z?-dimensional analogue of the usual vector H 
in 3-dimensional electromagnetism. As explained earlier 
it is necessary to treat the magnetic field as a tensor if 
D ^ 3. For more details of this analysis see reference 



We have treated the magnetic field gener ally in (|3lo|— 
3.43 ). Specializing now to a single component field of 
strength B {F12 — — -F21 = B), equation ( [3.43| ) can be 



written as 



where 



H = B- NL 



M 



Jb 



(3.44) 



(3.45) 



and H^"^ = —H^^ = H. Equation (3.44) is the conven- 
tional B-H relation found in three spatial dimensions, 
but with the notation defined here it can be seen to 
hold for all D. M in ( 3.45| ) is the magnetization. This 
approach is seen to avoid any ambiguity between what 
Schafroth called the actin g and microscopic fields. 

We can now split M in ( 3.45) in to two pieces in an ob- 
vious way. The derivative in ( 3.45 ) is a functional deriva- 
tive, and because B \s a. constant for our problem we can 
define 



1 917(0) 

V dB ' 
1 d^T^t) 

V dB ' 



(3.46) 
(3.47) 



Mt^o is easily computed using (3.24) to be 



TO \2'kI3 J 
-xS, [2-D,l]} 



(3.48) 
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Recalling (3.25) we have 



M' 



^x^,[2- D,l]\ . (3.49) 



In cases where no phase trans ition occurs, we have Qi = 
Q so that the first term of ( 3.49 ) is constant. Also if 
^ = the n M *^"^ — 0, so that the total magnetization is 
given by ( |3.49D with Qi^Q. 

When a phase transition does occur w e n eed M^^\ 
With our gauge choice ( 3.26| ), and using (U) for 
we find 



mV 



d"xy'^\^\ 



(3.50) 



was given by ( |3.9| ) with (3.27). It is easy to show that 

1 



which results in 



d"xy^\fo\' 



\Co\ 



2eB' 



2mV 



Qo 
2mV 



(3.51) 



(3.52) 



after using ( 3.1l| ) . We may now combine ( 3.49 ) and 
(|^ to read 



M 



Q 



e f m \ 
2mV ' m \2ttJ3 J 
-xE, [2-D,l]} 



(3.53) 



since Q = Qo + Qi is the total charge 



The result in (3.53) is the exact expression for the mag- 
netization which holds even if there is a phase transition. 
Had th e con densate ^ bee n ign ored, we would have ob- 
tained ( 3.49 ) rather than ( 3.53 ). In the true expression 
( |3.53D the first ter m is constant, whereas in ( |3.49| ) the 
first term is not constant if there is a phase transition. 
Thus, neglect of the condensate for I? > 5 would lead to 
an erroneous result for the magnetization. 

A selection of graphs for the dimensionless magneti- 
zation Ai = mXM. is shown in figure ^. Although the 
expression for the magnetization will not be correct for 
D > 5 if the condensate is neglected it can be seen that in 
the case D = 5, the graph is smooth and continuous for 
all values of the field. For D = 7 a kink is apparent. This 
seeming incongruity can be explained by the fact that 
the magnetization is calculated from a first derivative of 
the thermodynamic potential whereas the heat capacity 
comes from a second derivative. Hence the magnetiza- 
tion curve is smooth in those cases {D ~ 5,6) where the 
heat capacity is continuous. 

The zero-field spontaneous magnetization plotted on 
the graphs is of the Schafroth |9) form 



M 



Q 

2mV 



(3.54) 



and it can be seen that as the field B ^ 0, this limit is 
recovered. 






(3//?:, 



FIG. 8. The dimensionless magnetization A4 = 



for 



the charged hose gas in 3, 5, and 7 spatial dimensions. As 
the number of dimensions increases, the diamagnetism due 
to macroscopic occupation of the lowest Landau level in the 
low temperature (high ^) region becomes more pronounced. 
As can be seen from the graphs, the general behavior of the 
magnetized gas is similar in cases when a phase transition is 
absent {D = 3) or present {D = 5, 7). The even dimensional 
cases {D = 4, 6) have been omitted for brevity as they merely 
interpolate between the odd ones shown above. 
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IV. ANALYTICAL EXPANSIONS 

So far we have presented mainly the resu lts ob tained 
from a numerical evaluations of the sums ( |3.23 ), since 
it is not possible to obtain exact results. However due 
to the form of the exponential involving e^'^, in order 
to obtain reliable numerical results directly from ( 3.23| ) 
it is necessary to include an increasing number of terms 
as X decreases. In this section we will discuss a reliable 
method for obtaining approximate analytical expressions 
for various situations when e and x are small by finding 
asymptotic expansions for [a, (^l- 

The basic method we will use here involves the Mellin- 
Barnes contour integral representation for the exponen- 
tial function: 



(4.1) 



Here c is a constant with 5Re c > so that the contour lies 
to the right of the poles of the Gamma function. This is 
basically equivalent to the method used by Robinson |^] 
to obtain asymptotic expansions for the Bose-Einstein 
functions. It was used for the 3-dimcnsional magnetized 
gas by Daicic and Frankel jl^, and has been used to 
discuss Bose-Einstein condensation in a harmonic oscil- 
lator potential in reference There are various ways 
in which (4.1) can be used to obtain expansions over a 
range of discussed in references p3| , ^ . We will 

content ourselves with the simplest presentation here. 

The details of how ( [4.l[) may be used for the class of 
sums (3.23) are given in Appendix |A| for the situ ations of 
intere st to us. The basic method is to use (4.1) to con- 
vert (B.23) into a contour integral. The contour may be 
closed in the left hand side of the complex plane and the 
result evaluated by the residue theorem. The even and 
odd spatial dimensions differ somewhat in the pole struc- 
ture of the integrand. The net result is an asymptotic 
series for S]„ [a, 8\ which can be used to approximate the 
specific heat, magnetization and other thermodynamic 
quantities. 



A. Critical temperature 

For D > 5 the magnetized Bose gas is characterized 
by a well defined critical temperature Tc which satis- 
fies (3.35). It is not possible to evaluate Tc in closed 
form. However we can use our asymptotic expansion of 
Si [2 — D,0] to obtain an approximate result for weak 
magnetic fields. We have from Appendix |^ 



Q~eV 



{CH(f) + fc4^)+---} (4-2) 



if only the two leading terms are included. This assumes 
<C 1. With the free Bose gas critical temperature 
defined by 



= eV 



C4 



we find 



1- ^ 



(4.3) 



(4.4) 



Since Xc is assumed small we see that (3c 
to show from (4.4) that 



Tc 



Tn- 



1 C« 



' D-2 
. 2 



) eB 



D U § 



Pq. It is easy 



(4.5) 



to leading order in ^ . This shows that Tc ^ Tq a,s B ^ 
0. Furthermore, for a fixed charge density, the critical 
temperature is lower when a non zero magnetic field is 
present. This is consistent with our earlier numerical 
results. 

A cruder estimate of Tc was given in refe rence ]39| ] 
which had the same linear behavior as in ( |4.5| ) but with 
a different numerical factor in front of B. Our result (4.5) 
is a special case of the multi-component magnetic field 
presented in reference 

It is possible to improve on the linear approximation 
of {i.l ) by working consistently to higher order in the ex- 
pansions. It is necessary to deal with D = b,6 separately 
from D > 6 because of the order of the terms retained. 
For D = 5 we find 



Tc 
To 



For D 



1 



Xq 



1 Cn. 



1 c. (I) 

lOOCl (I) 30C„ 



6 we find 

^2 



57r2 Cr (I) " 

CH(i) 



o 



(4). (4.6) 



Tc TT" 

36G (3) 



Xo 



6C« (3) 



X 



C'(2) 



TT 



^2 216Ch (3) 
For D > 7 we have 



Xq 

In- 6 

2tt 



(4.7) 



1 Cu 



D-2 ^ 

2 J 



3(L» 



2D 



Xq + 



1 c. 

6^ C. (#) 



D-2 ^ 
2 / 



(f ) 



(4.8) 



It is worth remarking that the result given by May pl| ] 
is only correct if it is taken to linear order in B. 

We can also obtain an approximate expression for the 
charge in the condensate w hen T <TciorD>5. When 
T <Tc we have e = 0, so ( 3.25D gives us 



Qi=eV 



m 

2^ 



xSi [2-Li,0]| 



e=0- 



(4.9) 
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By using ( 3.35| ), which defines the critical temperature 
Tc7, we find (noting x = Puj, Xc = f3c^) 



D-2 



(4.10) 



From ( 3.10 ) wc find that the charge in the condensate is 



Qo = 



1 - 



[2-i^,0]|,^o 
[2- DM 



(4.11) 



This result is exact. If we now use the approximate ana- 
lytical expressions for Ei we obtain 



go = Q 



(a; - Xc) 



(4.12) 



It is worth remarking that the accuracy of any of our 
approximate expressions can be increased by simply in- 
cluding more terms. 

From ( 4.12| ) we see that as w ^ because ^ = 
and we know that Tc ^ To, we recover the free Bose gas 
result 



!o = 



1 - 



(4.13) 



The term in (4.12) which involves {x — Xc) represents the 
lowest order correction to the free field result in a weak 

magnetic field. 

In figure ^ we have plotted the free gas result (4.13), 
the exact numerical result, and our approxi matio n (4.12) 
for the case D = 5. It can be seen that ( 4.12 ) is very 
close to the true result. 




FIG. 9. The ratio of the charge containe d in the ground 
state over the t otal charge for i) the free gas ( 4.13 ), ii) the ex- 



act result (4.11), and iii) the approximation (4.12) (for D = 5, 

LO = 1). 



B. Magnetization for D > 5 



The magnetization may be evaluated using ( ^j.53| ). It 
is necessary to distinguish the cases D = 3,4 for which 
no phase transition occurs from D > 5. We will deal with 
D > 5 in this section, leaving D = 3, 4 until later. With 
D > 5 we must deal with T > Tc and T < Tc since the 
chemical potential is different in the two cases. A further 
complication is that even and odd dimensions must be 
treated separately, as explained in Appendix ^ and in 
addition the leading terms are different when Z? = 5, 6 
than for D > 6. For simplicity we will only give the 
expressions for T < Tc- When T > Tc we have e ^ 
and the results are more complicated. In addition, the 
expansions break dow n onc e e grows too large. 

For T <Tcwe use (|3.53D with e = 0: 



M 



'2mV m \2t:(3 J 
-xE, [2-D,l]}|,=o- 



{E, [~D,0] 



(4.14) 



This can be rewritten if we eliminate (^) ^ using (B.35) 
We find 



M 



Q 



2mV 



(El [-D,0]-x^, [^-D,l]) 



e=0 



Xc^. [2-D,0\\e=0. 

x=xc 



(4.15) 



It is now a straightforward matter to use the results of 
Appendix ^ to expand the sums. 



1. D = 5 



From (All) with a = —5 and e = we have 



2 12 



47r5 



720 



(4.16) 



From (A23) with a = —3, £ = and S = 1 we have 

Gil) Uil) SxiCnil) 



-3,1] 



12 



-'(^ (i) 

9607r 



(4.17) 



From (All) with a — —3 and e = we have 
^ ^« (f ) , ^« (f ) I ^cCj? (I) 



Si[-3,0]|. 



27r5 



c4i 

12 



18) 
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It is now easy to show that the asymptotic expansion of 
the magnetization is given by 



M 



where 



2mV 



1 



(4.19) 



27r2 G (2; 

^ ^-R (2) -^cC-R (2) 

1 6" 



^Xcxi Cr (I) 



CI (I) 12G (I) 

cl(i) 



+ 



By taking B 
magnetization 



87r2 Ch (f ) ^ ^ 

-> we are left with the spontaneous 



M{B ^ 0) ~ - 



2mV 



(4.21) 



since Tc Tq and f5{x,Xc) in this hmit. This 
is the 5-dimensional version of Schafroth's result, and 
agrees with May When B there are corrections 
to the Schafroth form as shown by the terms in /s (x,Xc)- 
There are two sources for these corrections. The first 
is that for B ^ 0, Tc ^ Tq. The second is that the 
asymptotic form of M has higher order terms present. 
The fact that M does not vanish as ^ shows that 
the Meissner-Ochsenfeld effect exists. 

We can also compare our result with tha t give n by May 
0. To do this we must replace Tc in ([4.19| ) with Tq. 
We will o nly work to first order in the magnetic field and 
use ( 4.19|) . It is easily shown that 



M ~ - 




(4.22) 



This agrees with May's result where a different method 
was used. There are two comments to make here. The 
first is that had we taken the expansion beyond linear 
order in x the results obtained would not agree with that 
of May's method for reasons already mentioned. Sec- 
ondly, care must be exercised in using (4.22) because we 
have assumed T < Tc in its derivation. This means that 
(4.22) does not hold a t T — Tq, as can be seen from [O. 



If we set T = To in ( |4.22| ) we would conclude that M 
was positive, leading to paramagnetic rather than dia- 
magnetic behavior. This is clearly wrong. We will return 
to results for T = To later. There is of course nothing 
wrong with taking T = Tc in ( l4.19| ). 



2. D = 6 



From 

S.[-4,0]U^o. 



(A9) with £ = 0, fc = 2 we have 
Ch(3) , Ch(2) 



C(2)^(ln|^-7) 



8640 



27r2 
+ 0(xl' 



12 



(4.23) 



From ([ASD with £ = 0, fc = 3 we find 

G(4) , C«(3) , xG(2) 



Si [-6,0] 



X 2 

^ x' 



12 



87r2 



1440 



(4.24) 



In this case all higher order terms in the asymptotic ex- 
pansion vanish. Finally from (A22) with e = 0, 6 — 1 
and k — 2 we find 

Again higher order terms i n th e asymptotic expansion 
vanish. In obtaining ( 4.23| — 4.25 ) various relations of the 
Riemann ^-function were used. The expression for the 
magnetization becomes 



where 

/e {x,xc) 



3 2 

3^2 



6C« (3) 



+ 



4^2 216CI(3) 144CI(3) 



G (3) 



C(2) 



+X-1 



TT 



27r2 

6 



In^ 
12 



1_ 
12 



6C2(3) 
2 



1728Ci(3) 

C(2) 
27r2 

6 



TT 



2592(3(3) 
C(2) 



18C2(3) 
^2 



In^ 
12 



In^ 



1_ 

12 



27r2 



12 



1_ 
12 



4Ch(3)V45 4 



(4.27) 



Again we see the Meissner-Ochsenfeld effect as B — > 0, 
with the magnetization of the Schafroth form. 
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3. D>7 



Using the results given in Appendix and working 
only to quadratic order in x, Xc (to avoid treating even 
and odd dimensions separately) we find 



M 



Q 



2mV 



where 



/d 



Since {x,Xc) 
magnetization. 



1- - 




Xc 

V3 2 



(f ) 



Cn 



' D-2\ 



a (f ) (f ) 
am , 



(4.29) 



6CI (f ) 

as _B ^ we recover the Schafroth 



3a: 



[(l-2e)vrCMi,^)CMi) 



-4<I(§)C.(i£)CH(i,e)]+0(xi). (4.30) 



2. L» = 4 

In even dimensions, as noted in the appendix, the anal- 
ysis is complicated by the presence of double poles. These 
give rise to terms in In x which dominate over the higher 
power terms which are the only ones present in the odd 
dimensional cases. In addition, the situation is compli- 
cated by terms involving Ine and i^{e). The expansion is 
given by 

C„ _ 36Ch (3) 27r2 (In a; - ipie)) 



l-2£ 



3(lna;)2 
216C„ (3) 



Ine + C;^ (0,1 + e) 



-O 



(4.31) 



C. Specific Heat 

In a similar manner to the previous section analytical 
expansions may be obtained for the specific heat capac- 
ities in the dimensions of interest to us in this paper. It 
is more convenient to express these results in the form 
they were given previously, i.e. as expansions of the ra- 



tio 



N 



These expansions are constructed from the ex- 



pression for the heat capacity (either equation ( p.3q ) if 
above the transition temperature, or ( 3.39| ) if below it) 
as well as that for the number density (obtained from 
( p5| ) as iV = f ) . The various sums present in these 
equations are treated as before by substituting the ap- 
propriate analytic approximations given in appendix W. 
Brief comments on the form of the expansions are made 
in each case. 



1. D = 3 

From figure || it can be seen that the charged Bose 
gas in three spatial dimensions does not exhibit a phase 
transition with finite field, but approaches the zero field 
result in the B ^ limit. The maximum of the heat 
capacity is always lower than the zero field limit, a fact 
which can clearly be seen from the analytic expansion 



Cy 



15Ch (I) 3x1 

4G(i) 40FcMi)c«(i,s) 

+5K«(l)CH(i,£)C«(i,^)] 



3. D = 5 

The five dimensional gas is the first to show the ex- 
istence of a phase transition and as has previously been 
discussed in subsection HI C , it is necessary to use dif- 
ferent expressions for the heat capacity above and below 
the critical temperature. For T > Tc, this quantity is 
given by 



a 35CH(i)G(l)-25CMI) 
^" 4C.(i)G(|) 

257T^xkH{le)Cn{l) 



^Cl (f ) 
5a;(2£- 1) 



sa (i) a (I) 

-2amim+o{xi), (4.32) 

while below Tc equation ( |3.39| ) must be used, and hence 
g ^ 35C. (i) X (l5Cg (I) - 35Ch (i)C. (i)) 

N 4c«(i) 8CMI) 

357r^:riC.(i)C.(-i) 
2CI (I) 

7x^(6C.(i)C.(i)+5C.(^)C.(i)) 



48CMI) 



+0{x^). 



(4.33) 
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I D = 6 

The D = 6 case is similar to the D = 4 case except 
that here too a phase transition occurs below a critical 
temperature. Again the presence of logarithmic terms 
complicates the expression. For T > Tc 

a 2 (tt^ - 405CI (3)) X 



N- 157r2G(3) 90^4^2(3) 
X [7r4(2e-l)(^«-135C^(3)) 
+29160CI(3) (C(0,l+e)-lne)] 
162(2£- l)a;lna; 

whilst for T <Tc 

a 27r4 X (270CI(3) - 7r6) n^x^ 



(4.34) 



iV ~ 15G (3) 90CI(3) 180CI(3) 

X [247r2C(-l,l)-2^2_^25G(3)] 

+0(x^lna;). (4.35) 



5. D^7 

D = 7 is the first case in which the heat capacity is 
wholly discontinuous. Above the critical temperature the 
expression 



77 



7(9C.(|)G(|) 



7CI(D) 



4CH(i)c.(D 

a:(14£ - 7) 



8c. (i)a (i) 



[5CI(DC.( 



+9G(|)c.(|)c.(l)-7cl(i){G(|) + c«(D}] 

©(x^), (4.36) 



49^5a;ic„(-i,e)C«(i] 



2C2 (f ) 
holds whereas below it 

C„ 63C.(|) , 7x(5Cl(i)-9C.(i)G(l)) 



TV 



8CI (i) 



[5C4i)c«(l)-35C«(l)c4i) 

(4.37) 



16CMI) 

-21CHa)CK(|)]+0(^') 



must be used. 



D. Expansions for T near To 

As noted in Schafroth's original paper ||^ the approxi- 
mation (for the 3-dimensional gas) 



M 



2mV 



(4.38) 



will break down when T becomes too close to Tq. This 
can be substantiated by a direct numerical evaluation 
of M and comparison with ( 4.3§| ) as we showed in p4| . 
The Schafroth criterion for validity of (4.38) can be de- 
rived in a simple manner as described in |17|. An obvious 
question to ask is that if ( |4.38 ) does not hold as T be- 
comes close to To, is there another simple approximation 
which can be used? This was first studied by Daicic and 
Frankel , who showed by expanding about ^ ~ that 
the magnetization M ~ —cB^ with c > a constant. 
This approximation was valid for values of T closer to Tq 
than Schafroth's approximation, but still broke down as 
T ^ Tq. In reference Q] we showed how it was possi- 
ble to evaluate the magnetization in a temperature range 
which included T = Tq. In this section we will present 
the details of this result and generalize the analysis to 
spatial dimensions D > 3. 

For D = 3 the critical temperature To for the free Bose 
gas is defined by 



Q^eV 



c.(i) 



where /3o = For B 7^ we have 

Equating these two expressions gives 



(4.39) 



(4.40) 



(4.41) 



where xq = f3oLO. The aim now is to solve this for s when 
the magnetic field is weak and T is close to To (meaning 
x is close to xq). Because we assume a weak magnetic 
field, we may use the first few terms in the asymptotic 



expansion of Si [ — 1, 



will approximate (4.41 



obtained from Appendix We 
) as 



C„(|) + (7ra;)kH(i 



e ] X ■ 



(4.42) 



Higher order terms could easily be included to improve 
the accuracy of the result. 

Suppose that we concentrate first on T = Tq. Call the 
value of £ at T = To, eo- Then (|0|) gives us 



o = CH(i 



£0 



£0 xk + 



(4.43) 



(The next term is of order x^). If we let i? ^ then 
Xq = Pqlo — > 0. Thus as ^ we must have £0 — > a 
where a is defined by 
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(4.44) 



The value of a can be found numerically with the result 
a = 0.302721829 .... We have verified this by solving 



(4.41) numerically for decreasing values of B and found 
that e — ^ a as S is reduced. Because /i = w — e) this 
result is still consistent with the expectation that /i — s- 
as i? 0. We have essentially determined how fast fi 
vanishes as B — > 0. 

For small, but non-zero, values of B we can try to 
solve (|4.43D . In order to obtain a consistent expansion 
from ( |4.43| ) it is fairly clear that we must have 



eo — a + aiXQ + a-zxo + 



(4.45) 



for some coefficients ai. 



0-21 



which can be found by 
substituting ( 4.45 ) into ( 4.43 ) and working to a consistent 
order in xq. It is easily shown that 



-aixgCn (|,a) +O{xo). 



Use of ( |I45| ) and ( |44^ ) in ( pl3| ) fixes 



ai 



2C. (i) 



2-" 



(4.46) 



(4.47) 



It should be clear how we can obtain an approximation 
for So to any order in xq by extending the procedure we 
have just described to higher order. 

So far we have just concentrated on the evaluation of 
e at the single temperature Tq. Suppose that wc now 
extend this to temperatures which are close to Tq. By 
a simple extension of the analysis just presented it is 
possible to show 



Consistency of this expansion requires 



(4.48) 



(4.49) 



In particular, the approximation is good sd x — xq- 

The approximation we have found for e may now be 
used in the expressions for the specific heat and the mag- 
netization. For the specific heat at T = To we find 



'n 



9G 



4^ (i) 



2vr5C«(|,a) 




9G(^)C«( 



27CH(^)CH(|)a(i,a) fl 



2<l {\ 



-O 



(4) 



(4.50) 



This shows that as B ^ the specific heat approaches 
the free field result at T = Tq, confirming analytically the 
trend we found numerically. It also provides an analytic 
proof that for small xt^, the presence of the magnetic 
field lowers the value of the specific heat from the free 
field value. 

For the magnetization we find (at T = Tq) 




(4.51) 



This provides confirmation of the M ~ —cB^ magne- 
tization law of Daicic and Frankel |l^, but at a lower 
temperature (and therefore the coefficient oi B^ differs 
from that of ||l^). In addition, we have computed the 
next order correction to the leading B^ behavior. 

Having presented the case D = 3 in some detail, we 
will now examine other dimensions in less detail. For 
D = 4, in place of ( 4.41 ) we have 



xSi [-2,0] = 
with Tq now defined by 

Q = eV 



Xo 



Cu (2), 



Cn (2). 



(4.52) 



(4.53) 



(Here xq = Pquj with = Tq ^ as before, and we note 
(2) = ^). Using the expansion for Si [—2,0] of Ap- 
pendix valid for small x and small ex we obtain 



: 



G (2) 



X In (ex) 



ex - 



(4.54) 



Concentrating on the temperature T = Tq at which e 
Eq we find to lowest order, that eo must satisfy 



= ( £o - ^ ) In {eoxo) 



£o- 



(4.55) 



Assuming Eq is finite as xq 
use I In I ^ I In eo I to find 



0, then in (4.55) we may 



1 



2 i ' Inxo 



(4.56) 



if only the leading term is kept. This shows that eo — > ^ 
from below as xq 0. Unlike the case D = 3, for 



which lim^^o 
find hms^o 4 



B 



= ^ — a was non-zero, this time we 
0. Again by including more terms it is 
possible to render the approximation for eo more accu- 
rate. Because of the presence of the logarithm in (4.56), 
the convergence of eo towards the value of ^ is very slow. 
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If we expand the magnetization for T = Tq using 
it is found that 



M : 



Q xo^^f^ 

xo 



2mV 7r2 



(4.57) 



for small xq (Only the leading term has been included 
here). For the specific heat we have 



a G(3) 4G(2) 
N - G(2)^ln(f) ^ 



(4.58) 



if only the leading order term is kept. Again we see the 
approach of the specific heat towards the free field value 
at T = To as i? ^ 0, and that the specific heat is lower 
when B 0. The analytical results are consistent with 
the numerical ones presented earlier. 
For D > 5 we may use 

xo^^[2~D,0]U„^CI^{i) (4.59) 

to solve for Eq- (We still use xq = Pqcu and Eq to denote 
the value of e at T = Tq). From our earlier discussion 
we know that Tq > Tc so that Sq ^ 0. The expansion 
for Si [2 — D, 0] dep ends upon whether D is even or odd. 
For odd D we use ( Al6| ) with a = 2 — D and i5 = to 
obtain 



= C« (^) [I - so) + (^) (el - so 



4-D 



C,{*^,So)xo^ +0{xl). (4.60) 



By letting xq ^ (corresponding to i? — > 0) we see that 



In order to see how fast eo ^ ^ as — > we 
must distinguish between D = 5 and D > 5. For _D = 5 

the next to leading order term is of order x§ and we find 



£0 



1 

71 



+ O (xo) . 



(4.61) 



For D > 5 the next to leading order term is of order xq 
and we find 



IMS 

> 2 / 



24 C« 



(4.62) 



(The next term in (4.62) is of order Xq if Z) = 7 and 
for D > 7). 

In the even dimensio nal c ase we may use ( |A15| ) with 



k = ^ and 5 = in (k^). For I? = 6 we find 



1 Xo In Xq 



47r2 



Oixo). 



(4.63) 



For D > 6 the result in (4.62) holds. Therefore for all 
£) > 5 we find that Eq ^ as B 0, proving that 5 — > 
as i? ^ 0. This contrasts for the result for D = 3. 



V. DISCUSSION AND CONCLUSIONS 

This paper has studied the thermodynamic properties 
of the ideal charged Bose gas in some detail. Spatial 
dimensions with D > 3 have been examined since even 
for the free gas it is known that the properties of the 
gas are sensitive to the spatial dimension. The specific 
heat was calculated numerically as well as analytically, 
for small values of the magnetic field. We also performed 
calculations of the magnetization and showed how the 
effective action method could be used to account for the 
condensate when D > 5. 

One motivation for our study was to understand in 
more detail the behavior of the magnetized gas for D = 3. 
When a magnetic field is present, no matter how small, 
there is no phase transition; however when there is no 
magnetic field the system does exhibit a phase transition 
with a non-zero condensate. On physical grounds we 
would expect that there should be some sense in which 
the system in a small magnetic field should behave in al- 
most the same way as the free gas. By examining the spe- 
cific heat it is possible to see that as the magnetic field is 
reduced, the curves start to resemble the specific heat for 
the free gas. So long as the magnetic field remains non- 
zero the specific heat is always smooth, with the specific 
heat maximum approaching the the free Bose gas tran- 
sition temperature as the magnetic field is reduced. The 
results of section IV may be used to study this analyti- 



cally. Although we did not show it, it is straightforward 
to show that the derivative of the specific heat becomes 
discontinuous as i? — > 0, exactly as in the case for the 
free Bose gas. Similar remarks apply to gases in other 
dimensions. 

We also studied the behavior in large magnetic fields. 
Here we found that the specific heat for a gas in D spatial 
dimensions looked like the specific heat for the free gas 
in (D — 2) spatial dimensions over a range of tempera- 
tures. Once the temperature gets too large this effective 
reduction in dimension disappears. 

By using the Mellin-Barnes integral transform we were 
able to obtain a number of analytical approximations. 
Although, for £> > 5, it is not possible to solve for the 
critical temperature exactly, it is possible to obtain good 
estimates when the magnetic field is weak. These ap- 
proximations can then be used to study the Meissner- 
Ochsenfeld effect. We were also able to obtain reliable 
approximations for the first time valid at the free gas 
condensation temperature. As a by-product our calcula- 
tions showed exactly how the chemical potential vanishes 
as B ^ 0. 

The most obvious way that the calculations given 
in this paper should be extended is by the inclusion 
of Coulomb interactions among the charged particles. 
Clearly this is essential before it will be possible to study 
Bose-Einstein condensation of charged particles in a re- 
liable way. Given the recent experimental advances in 
the cooling and trapping of atomic gases, it may become 
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increasingly important to study this problem for trapped 
ions. A less pressing extension of our work is to relativis- 
tic charged particles. Daicic and Frankel have 
done a study of this already in various cases; however, 
it would be easily possible to extend the analysis of our 
paper to obtain the specific heat for the first time. 
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APPENDIX A: DETAILS OF THE 
MELLIN-BARNES TRANSFORMATIONS 

1. K = 1 case 



The sum Ei [a, 6] defined in (3.23) is 

g-Za:(e+(5) 

1=1 ^ ' 



(Al) 



By using the binomial expansion on the denominator we 
find 



oo oo 



El [a, <5] = ^ ^ g-ix(n+e+5) 



(A2) 



1=1 n=0 



This is now in a form where (4.1) can be used for the 
exponential. There are a variety of ways to use (4.1) on 
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( [A2| ) , depending on how we split up the exponential. We 
require results for S = 0,1. 

Suppose that we concentrate on ^ = initially. If we 
are interested in the £ — *■ limit, it p roves advanta geou s 
to separate off the n — term in ( A2 ) before using ( [4.l[) . 
This is a slight variation on the method of reference ||17|] 
which is useful for isolating terms which diverge as e — )■ 0, 
although we do not have to do this. We find 

E,[a,0]-Li_„ [e— ]+ / :^r(%~« 

Jc — ioo 

xC«(e-f)C«(^,£ + l) (A3) 

if it is noted that the sums over I and n can be done in 
terms of the Riemann and Hurwitz C-functions respec- 
tively. If the n = term is not removed, then in place of 
(A3) we find just the integral part with e) in place 

of (h {S, 1 + e). We take c > real and large enough so 
that the contour in ( |A3| ) lies to the right of the largest 
pole of the integrand (specifically, c > max (1,1 + f), 
where max denotes the maximum of the set). The con- 



tour in (A3) is closed in the left hand side of the com- 
plex plane, and the residue theorem used to generate the 
asymptotic expansion of Si [a, 0] from the known prop- 
erties of the F- and ^-functions |Q . 

At this stage it is necessary to identify a number of spe- 
cial cases. If a = 0, then the poles of the Riemann and 
Hurwitz (^-functions both occur at = 1 and the inte- 
grand has a double pole. If a = —2k where fc = 1, 2, 3, . . ., 
then the pole of the Riemann ^-function coincides with 
one of the poles of the F-function and again the integrand 
has a double pole. For a ^ 0,-2,-4,... the functions 
occurring in (A3) all have distinct poles, so that the in- 
tegrand only has simple poles. 



a. Q = 0, 5 = 



Setting a = in (A3) results in Lip [e ^^]. This is 
simply 



Lio [e"^1 



(e-^^ - 1) 



(A4) 



which follows directly from the definition (^^). Evalua- 
tion of the integral by residues results in 



Si [0,0] 



1 



+ e) + Ina; 



(e^^ - 1) 



xC«(-l-2p,l + £). 



(A5) 



We use the symbol ~ here to denote that the result is 
an asymptotic expansion. The Riemann and Hurwitz 
functions in (A5) may be related to the Bernoulli num- 
bers and polynomials Ea]. The divergence as e is 



contained in the first term. When ex is small, we may 
use ra 



1 1 1 ^ (£x)^P+i 

(e^^ -1) ^ ex 2 ^ F(2p + 2) 



C„(-l-2p) (A6) 



p=0 



to obtain 



Si [0, 0] ~ — - i - i (VXl +£) + Inx) 
ex 2 X 



ix)^P+^ 

b;r(2p + 2) 



-E 

p=0 



G (-1 - 2p) 



{e'P+'+CH{-l-2p,l + e)}. (A7) 



b. a ■ 



-2k, fc = 1,2,3, . . .; <5 = 



We set a — —2k in (A3) and note that the integrand 
has simple poles at ^? = 1 and 6 = —p for p 7^ fc — 1, and a 
double pole at = 1 — fc. It is straightforward to expand 
about these poles and use the residue theorem to obtain 



Si [-2fc,0]~Li, [e--^]^ 
(-a;)'=-i 



Cn (k + 1) 



[C{l-k,l + e) 



m 

+ (7 + tp{k) - lnx)Cff (1 -k,l + e) 



p=0, 
pjik-1 



r(p + i) 



xCh (-P, 1 + e). 



(A8) 



Here (^'^{1 — k, 1 + e) denotes the derivative of {s, 1 + e) 
with respect to s evaluated at s = 1 — fc. For fc = 
1, Cf/(0, 1 + e) niay be found in It is possible 

to generalize the procedure described in |Q to obtain 
C^(l — fc, l-|-e) for other values of fc , although we have not 
found it possible to obtain as simple a result for fc = 1. 
The values we need are given in Appendix 

The polylogarithm Ijifc [e"*^^] may be expanded for 
small ex, also using (4.1). This was originally described 
by Robinson Q. We find 



Li. 



(-ea;)*^i (i/)(fc) + 7 - In (ex)) 



F(fc) 



E 

/=0; 



(-ex)' 
T{1 + 1) 



Uk-l). 



(A9) 



Finally we note that 



V'(fc) -f 7 = 1 + ^ + • • • + 
for fc > 1, and vanishes for fc = 1. 



(AlO) 
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c. q/0, -2, -4, ...;(5 = 

In this case all of the poles of the integrand in ( |A3| ) are 
simple, and we obtain 



Si[a,0] ~Li_. [£ 

+rfi 



C«(i-f) 



|)x~(i+t) 
{-xY 



C«(l + -,l+e 



xCff +£)• 



(All) 



The expansion of the polylogarithm for small ex is again 
described in reference E2] and is 



Li„^ [e- 



'^]^r(i + -)(ex)-i-t 



(A12) 



We also require the expansion of [ct,S] for S > 0. The 

difference between 6 = and (5 7^ can be seen from 
(A2): when 6 = the n = term can lead to a divergent 
sum over I when e — > 0. For 6 > the sum over I is 
always convergent, even when e = 0. Thus for 6 > Owe 
do not h ave to separate off the n = term . Using (4.1) 
for (A2) we obtain 



-^r{9)x-\,{e-^)(:^{0,e + 6). 

(A13) 

The same cases arise as for 6 = 0, and the analysis is 
similar. We will simply list the results. 

d. a = 0, S> 



Si [0,(5]- iyj{e + 6)+lnx) 

X 



E 



X 



2p+l 



r(2p + 2 



El [~2k,6' 



p=Q 

xC«(-l-2p,£ + <5). 



e. a = -2k, k = 1,2,3, S > 
Cn {k + 1) 



(A14) 



{-xr 



[C{l-k,e + 6) 



m 

+ (7 + ^(fc) -lna;)C„(l- fc,e + (5)] 



p=0, 



r(p + i) 



El [a, (5] 



xC„ {-p,e + 6). 
/. a / 0, -2, -4, . . . 

c«(i-f) 



(A15) 



5 > 



.r(i.|).-(-t)c.(t.|,.- 



+ l^„r(p + i)^n 

x(H{-p,e + 6). 



(A16) 



One final comment we wish to make is that ( [A14| — 



AI6 ) hold even for 6 = 0, although as we explained, in 
this case it may be more difficult to isolate the terms 
which diverge as e ^ 0. Nevertheless, these results can 
be used if desired. 



2. K = 2 case 



The sum E2 [a, 6] defined in (3.23) is 



E.[a,5] = ^- 
'-1 U 



;f g-/x(e+i5) 



(A17) 



The analysis of this sum proceeds in the same way as 
for El [a, (5]. Binomial expansion of the denominator of 
( |A17] ) gives 



E, [a,<5] = ^^;t(n + l)e 



— lx{n-\-e-\-5) 



(AI8) 



/=1 n=0 



We only require the result for 6 = 1 in the main part of 
the paper, so that the sum over I is convergent even for 
£ = 0. We do not separate off any terms. Application of 
(O) to (|A1^) resuhs in 



E2 [a, 6] 



:r(0)x-%(e-|) 



-too 27rr 
+ (l-£-<5)C„ {e,e + 5)}. 



(A19) 



(The sum over I is a Riemann ^-function, and the sum 
over n can be done in terms of the Hurwitz ^-function 
as shown). 
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a. a = 



The integrand of ( A19D has a simple pole at 9 ~ 2, 
a double pole at 6* — 1, and simple poles at 6 = 
0,-1,-2,.... We find 



4 

"1 
2 



2' ' 24 
-e-5-(l-£ - (5)V'(£ + (5) 



-(l-£-(5) Inx] 

x{C«(-2-2p,e + (5) 

+ (l-e-5)G(-l-2p,£ + ,5)}. (A20) 



r(p + i) 



p=0, 

X {C^j l,£ + 5) 

+ (l-£-J)C^(-p,£ + <5)}. (A22) 



d. Q 7^ 2,0,-2,-4,... 

The integrand only has simple poles in this case, and 
we find that 



s.[a,^] = r(i + |).-(^+f){c„(f,£ + ^) 

+ (l-£-<5)C.(l + f,e + <5)} 
C«(2-f) , (l-e-5XH(l-f) 



6. a = 2 



The integrand of (A19) has a double pole at = 2, 
and single poles at0 = l,O,— 1,.... The residue theorem 
gives, 



[2, ,5] 



2a; 

+ [1 - Inx - i/'(£ + (5) 

+ (l-£-^)C„(2,£ + 5)] 

x{C„(-l-2p,£ + (5) 

+ (l-£-J)C^(-2p,£ + <5)}. 



c. Q = -2fc, fc = 1,2,3, 



(A21) 



The integrand of ( A19 ) has simple poles at 9 — 2,1, 
a double pole at 6* = 1 — fc, and simple pole at 6* = — p, 
p = 0, 1, 2, 3, . . . with p 7^ fc - 1. We find 



Ch(2 + A;) , (l-£-(5)G(l + fc) 



-(7 + V'(fc) ~ In a;) 



r(fc) 

X {C„(-fc,£ + 5) 

+ (l-£-5)C^(l-fc,£ + (5)} 

+ (l-£-5)C(l-fc,£ + '5)} 



■E 



r(p + i) 



c.(-.-f) 



p=0 

X {C„(-p- l,£ + (5) 

+ (l-£-5)C^(-p,£ + ^)}. (A23) 



e. AC = 3 case 



From (3.23) we have 



S3 K'Jl^E 



g-/x(e+5) 



(A24) 



Again we are only interested in (5 > 0. Following the now 
familia r st eps of binomially expanding the denominator, 
using (4.1) for the exponential, and performing the sums 
in terms of C~functions results in 



E3 [a, S\ 



27ri ^ 2 ^ V 2 / 



x{C„ (0-2,£ + J) 

+ (3 - 2£ - 2(5)C„ (61-1, £ + 5) 

+ (l-£-(5)(2-£-<5)C«(0,£ + 5)}. (A25) 



/. Q = 



The integrand in ( A25 ) has simple poles at 9 
3, 2, 0, —1, —2, . . . and a double pole at = 1. We find 



C«(3) , (3-2£-2^)C«(2) 



1 

2x 



2a;2 
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-e-S)i2~e 
{CH{-2,e + S) 



S){^{e + S)+lnx}] 



i. a ^4: -2k, fc = 0, 1,2, . 



+ il-e-S){2-e-S) Q - £ - 

X {C„ (-3 -2p,e + S) 
+ (3 - 2e - 2S)C„ (-2 -2p,e + 5) 
+ {l ~ e - S){2 - e - S) 
xC«(-l-2p,e + J)}. 



g. a = 2, A 



(A26) 



The integrand of (A25) has simple poles aX 6 = 
3, 1, 0, —1, —2, . . . and a double pole a± 9 — 2 for a — 2. 
When a = 4 we have a double pole at = 3, and simple 
poles at = 2, 1, 0, . . .. Because we do not require these 
cases in our paper we will not give the results here. 



In this case all the poles of (A25) are simple and we 
obtain 



+ (3-2e-2(5)C„(|,e + 

+ (l-e-J)(2-e-^)C„(l + |,e + <5)} 
CH(3-f) (3-2e-2<5)C«(2-f) 



+ 



+ 



2x2 



{l-s-5)(2-e-5)U{l-%) 



2x 



E 



[-xY 



2Y{v 



1)^4-^+2) 

p=0 

x{C„(-p-2,e + <5) 

+ (3-2e-2(5)C„(-]3-l,e + ^) 

+ (1 - e - ^)(2 - e - ,5)Ci^ (-p, e + (5)} . (A28) 



APPENDIX B: DERIVATIVES OF THE 
HURWITZ C-FUNCTION 



The integrand of (A25) has simple poles at 6* = 
3, 2, 1, —p for p = 0, 1, 2, . . . but p ^ k — 1, and a double 
pole a,t 6 — 1 ~ k. We find 



'2k, 5] 



Ch(3 + k) 
x^ 

(3-2£-2(5)C«(2 + fc) 



2a;2 

(l-£-<5)(2- 



2x 



{-xf 



fc-i 



{[7 + ^(fc)] 



2r(fc) 

+ (3-2£-25)C«(-fc,£ + 5) 
+ (1 - £ - (5)(2 - e - 5)C« (1 - fc, £ + 5)] 
+C(-l-fc,£ + 5) 
+ (3 - 2£ - 25)C (-fc, £ + 5) - In a; 
+ (1 - £ - (5)(2 - £ - 5)C (1 - fc, £ + 5)} 
(-xf 



E 

p=0. 



2T{p + 1 



■G(fc-P) 



The derivative of (s, a) with respect to s at s = is 
given by El 



C(0,a)-ln 



Via] 



(Bl) 



The derivative makes use of the Plana summation for- 
mula to derive an integral representation for Qh (s, a) (the 
Hermite representation) which is then differentiated with 
respect to s. In a similar manner the derivatives at other 
values of s may be calculated, although the results are 
less simple than (IbTI). For s = — 1 we find 



(6a2(7 + 1) + 6a + (1 - 7 - In 2ti)) 



12 



In 



na)V\ , C(2) 



2tt } 



+ 



27r2 



t=i 



it + 2) 



Here Cr(2) is the derivative of the Riemann (^-function 
given by 



C(2)--E 



Ini 



(B3) 



n=l 



x{a(-p-2,£ + 5) 

+(3-2£-2(5)C«(-p-l,£ + J) 

+ (1 — £ — (5)(2 — £ — (5)Cff (— p, £ + (5)} . (A27) which can be evaluated easily numerically. For s — —2 

we find 
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, ^_ (a3(87 + 6) + 9a^ + a(3-27-21n27r)) 

+ini yiMrV^^"^^^^^"*^"^^^ 



t=i 



2^1 y V 47r2 

'a(*+3)ft + 4)C 
(i + 2)(i + 3) 



°2, (_l)(«+i)a(t+3)(^ + 4)^^(^+l) 

-I- 9Vf -I- .-^^ ■ ' 



A useful check on the results is provided by the general 
relation 

— a) = -C„ {s + 1, a) - sCis + 1, a). (B5) 
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